High-energy antiprotons from old supernova remnants 
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A recently proposed model [lj explains the rise in energy of the positron fraction measured by 
the PAMELA satellite in terms of hadronic production of positrons in aged supernova remnants, 
and acceleration therein. Here we present a preliminary calculation of the anti-proton flux produced 
by the same mechanism. While the model is consistent with present data, a rise of the antiproton 
to proton ratio is predicted at high energy, which strikingly distinguishes this scenario from other 
astrophysical explanations of the positron fraction (like pulsars). We briefly discuss important 
implications for Dark Matter searches via antimatter. 



PACS numbers: 98.70.Sa 

Introduction — The antimatter component in cosmic 
rays (CRs) has been recognized since long time as an 
important diagnostic tool for cosmology (e.g. matter- 
antimatter asymmetry of the local universe), particle 
physics (indirect dark matter searches), and properties 
of cosmic ray sources and propagation medium (see the 
standard textbooks @ or the more recent reviews [1,0]). 
Recently, the PAMELA satellite detector Q has pre- 
sented its first results of the measurement of the positron 
fraction in the cosmic ray spectrum, which appears to 
begin climbing quite rapidly between ~ 7 GeV and 100 
GeV @ . This trend confirms (with much higher statistics 
and over a wider energy range) what previously found by 
other experiments, including HEAT Q and AMS-01 [If. 
On very general grounds, this behavior is at odds with the 
standard predictions for secondary positrons produced in 
the collisions of cosmic ray nuclides with the inter-stellar 
medium (ISM); an additional source of positrons seems to 
be required [9( . While numerous models of dark-matter 
(DM) annihilation or decay have been proposed (for a 
complete list see refs. to Q), astrophysical explanations 
do exist, in particular invoking e + — e~ acceleration in 
pulsars flO ]. 

A generic feature of these astrophysical solutions is the 
absence of a significant anti-proton signal accompanying 
the positron one, since the acceleration involves purely 
electromagnetic phenomena in pulsar magnetospheres. 
Actually, the PAMELA collaboration has presented also 
data on the p/p ratio below E ~ 100 GeV [HI, which fit 
naturally in a scenario of purely secondary production 
via CR spallation in the interstellar medium (ISM). In 
turn, this puts a non-trivial constraint on dark matter 
models trying to account for the positron excess [l2| ] . 

Recently, one of us has proposed an alternative and 
even simpler astrophysical explanation for the feature 
observed in the positron fraction PJ. In this scenario, 
the 'excess' is due to positrons created as secondary 
products of hadronic interactions inside the standard 
sources of CRs, supernova remnants (SNRs). In particu- 
lar, positrons would be produced in the late stage of SNR 
evolution, when also the bulk of cosmic rays (namely be- 
low the knee) are expected to be accelerated. The crucial 



physical ingredient which leads to a natural explanation 
of the positron flux is the fact that the secondary pro- 
duction takes place in the same region where cosmic rays 
are being accelerated; secondary e + (and e — ) participate 
in the acceleration process and turn out to have a very 
flat spectrum at high energy, which is responsible, af- 
ter propagation in the Galaxy, for the observed positron 
"excess" . The values of the parameters which lead to an 
explanation of the rising positron fraction are typical of 
old SNRs, rather than the young, often gamma-ray and 
X-ray bright ones. Since this is now a hadronic mech- 
anism for the explanation of the data, one expects an 
associated feature in the antiproton spectrum. The pur- 
pose of this letter is to present a preliminary calculation 
of the p/p ratio within the simple model of [l| . It is im- 
portant to realize that this model applies to a stage of 
the SNR evolution in which: 1) not many observations 
are available, with the possible exception of the ones in 
the radio band, 2) many effects are expected to play a 
role, such as magnetic field damping, on which we have 
exceedingly poor control, and 3) it would be important 
to carry out the calculations in a time dependent way, 
in order to move beyond a simple estimate. Thus, some 
of the parameters adopted in [l| might be considered as 
"effective" astrophysical inputs. While they need to be 
checked versus more realistic models, for the time being 
we believe that a more urgent task is to establish whether 
the model of [![ can account for the rising trend of the 
PAMELA data; this can be done most reliably by check- 
ing it versus independent predictions of the same model. 
The calculation of the antiproton flux (or, rather, the an- 
tiproton to proton ratio) is the first of them and the one 
requiring the minimum number of independent assump- 
tions. Such a calculation reveals that: i) the additional 
signal does not violate existing data; ii) a generic predic- 
tion of the model is a flattening and eventually a weak 
rise of the p/p ratio in the decade ~ 100 ~ 1000 GeV; 
the value of this ratio at TeV energy is about one order 
of magnitude above expectations from the conventional 
models (of course, this is equivalent to saying that a spec- 
tral break is predicted in the absolute p spectrum). Since 
this feature is strictly related to the one in the positron 
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spectrum, the model is very predictive. Clearly, these 
results have very important implications also for dark 
matter searches via anti-protons as well as for astrophys- 
ical diagnostics via the p/p ratio, as we shall comment at 
the end of the letter. 

The calculation — Here we only report the specific equa- 
tions that are needed to calculate the p/p spectrum, while 
referring to [l| for a detailed description of the model. It 
is also worth stressing that up to several tens of GeV 
the p/p spectrum is well in agreement with conventional 
mechanisms [TTj . Since we are only interested in the high 
energy part of the spectrum, E ^> 10 GeV, we shall im- 
plicitly assume approximations valid in the relativistic 
limit. Compared with the treatment for positron produc- 
tion proposed in ref. [l[, the differences arise only in the 
production cross-section and in the physical processes rel- 
evant for propagation. Concerning propagation, for the 
purposes of this letter we can neglect energy losses, re- 
acceleration, tertiary production, solar modulation, etc. 
which are relevant at low energy. The antiprotons are 
injected inside the sources as described by the function 

Qp(E) ~ 2 / °* X d£N CR (£)a pp (£, E)n gas c , (1) 

JE 

where c is the speed of light, n gas is the gas density for pp 
scattering in the shock region and a p p(£, E) is the differ- 
ential cross section for a proton of energy £ to produce a p 
of energy E. The energy £max is the maximum energy of 
the protons being accelerated in the SNR at the age rele- 
vant for the mechanism discussed here and it is discussed 
in the following. The factor 2 accounts for the antiproton 
coming from antineutron production, which we assume to 
be identical to the p one (isospin symmetry limit). A sub- 
tle point is that antineutrons, being neutral, stream freely 
away from the acceleration region until they decay (bar- 
ring nuclear collisions). The range of an (anti-)neutron of 
energy E is R n ~ (E/m n c 2 )r n c ~ 10~ 5 -EgcV pc, where 
T n and m n are the lifetime and mass of the neutron, re- 
spectively The following considerations assume that the 
confinement/acceleration region has a characteristic size 
I 3> R n , which is easily fulfilled for pc-scale shocks in 
SNRs (see also below). As in most calculations in the 
modern literature, for the cross section <jp(£,E) we use 
the parameterization of Ref. After production, the 
spectrum described by Eq. (fT]) is modified by accelera- 
tion inside the source and by propagation to the Earth. 
The latter phase is identical for p and p, so the spec- 
tral modification induced by propagation cancels in the 
antiproton to proton ratio. We assume that SNRs ac- 
count for the overall CR flux at the Earth and, for the 
moment, we are assuming that it is entirely made of pro- 
tons. Also, throughout the paper we are relying on the 
fact that most of the GeV-TeV production of cosmic ray 
protons happens in the late stage of SNRs which is of 
concern here. Then, the solution for the p/p flux ratio 
can be easily derived from [l[ in the form 

SNRs(E) 



where 



A(E)=-y(j + f*) x (3) 
x /"^dw^- 3 ^ d££ 2 ^a pP {£,uj),(4) 
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d££ 2 ^ a pp {£,E). (5) 



In the above expressions, n\ is the background gas target 
in the upstream region of the shock, u\ the fluid velocity 
there (which we fix at u\ = 0.5 x 10 8 cm/s), t$n is a 
typical SNR age, here fixed to t$n = 2 x 10 4 yr. The 
parameter £ (which we fix as £ ~ 0.17) is the fraction of 
proton energy carried away by a secondary antiproton, 
while r is the compression factor between upstream and 
downstream. The index —7 is the slope of the spectrum 
in momentum space, which is related to the spectral index 
a in energy space of the accelerated cosmic ray protons 
at the source via a = 2 — 7 and to the ratio r by 7 = 
3 r/(r — 1). Here, we fix r = 3.8 so that a ~ —2.07. Note 
that this is another instance of the oversimplification we 
are forced to here: the compression factor r is chosen in 
order to achieve an injection spectrum oc E~ 21 necessary 
to fit the CR spectrum after propagation. However, it 
is well known that SNR shocks stay strong (r = 4) at 
almost all times. It follows that the spectrum steeper 
than E~ 2 should follow from a complex overlap over time 
during the SNR evolution rather than the fact that the 
shock is weaker. If taken into account, this effect leads 
to an enhanced rate of production of secondaries inside 
the SNR for a given set of parameters. 

Finally, the function D\(u>) is the diffusion coefficient 
upstream of the shock, which in quasi- linear theory writes 
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where B is the magnetic field, e the unit charge, T ~ 
{AB/B) 2 is the ratio of power in turbulent magnetic 
field over that in the ordered one, A c is the largest co- 
herence scale of the turbulent component, and f3 is the 
index carachterizing the fluctuation spectrum. For sim- 
plicity, in the following we fix (3 = 1 (Bohm-like diffusion 
index), in which case the model does not depend explic- 
itly on A c (for a Kraichnan spectrum, one would have 
f3 = 1/2 and a dependence from the square root of both 
E and A c ). Denoting by B p g the magnetic field in micro- 
Gauss and by EceV the energy in GeV, numerically one 
has 



Di(E) 



3.3 x lO 22 ?- 1 E Ge vB 



uG cm2 s ~ 



(J) 



J p (E) 



2n 1 c[A(E)+B(E)\ 



(2) 



For the following numerical estimate, we fix n\ = 2 
(in cm -3 ) and {J 7 , B p q} = {1/20, 1}. We consider these 
numbers as reasonable if applied to old SNRs, in which 
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magnetic field amplification is not effective and in fact 
it is likely that magnetic fields are damped (see for in- 
stance [17j), We stress once again that this period is 
very poorly modeled and a precise quantification of the 
astrophysical parameters is tricky: for instance damping 
is required to lower the maximum energy of accelerated 
particles, but the temporal dependence of the maximum 
energy is not known, though it is expected to be rather 
fast. The velocity of the shock u\ is better known, since 
it can be estimated by using the standard Sedov solu- 
tion in a constant density of the background medium, 
yet the new term is quite sensitive to it (depending on 
u\ ). More complicated situations — such as the expansion 
in a density profile induced by a presupernova wind — are 
of relevance only in the early stages of the expansion of 
the shell, and in any case only for supernovae of type II. 
All in all, we are using simple effective parameters with 
all the limitations that this approach implies. More im- 
portant for the phenomenology is that the combination of 
parameters n\B~QU^ 2 / T « 160 is roughly what required 
to fit the high-energy behavior of the positron fraction, 
within a fudge factor of 0(1). Note that, for the chosen 
parameters and the energy range we are interested in, the 
characteristic size of acceleration I ~ D\/u\ is roughly 
three orders of magnitude larger than R n , confirming a 
posteriori the validity of including n in the source term. 
Another important point to discuss is that of the max- 
imum energy for the primary and secondary particles: 
protons accelerated at the shock have a maximum en- 
ergy which in principle can be estimated by equating 
the acceleration time and the age of the remnant. The 
maximum energy of secondary products is determined by 
the process responsible for their production: for electrons 
and positrons, typically the energy of the secondaries is 
£ ~ 0.05 of the parent proton. For antiprotons this frac- 
tion is £ ~ 0.17. However those secondary particles which 
are produced within a distance of order D(E)/u on both 
sides of the shock participate in the acceleration process 
and they end up being accelerated at roughly the same 
maximum energy as the parent protons, with a rather flat 
spectrum. For typical values of parameters one can easily 
find maximum energy in the range between 3 TeV (for 
Bohm) and 80 TeV (for Kraichnan). However these num- 
bers do not take into account a number of phenomena, 
such as the presence of higher energy particles generated 
at previous times, damping of the field and the possibility 
of obliquity of the magnetic field lines over most of the 
shock surface. Because of these numerous uncertainties, 
we adopt a sort of effective value of 10 TeV for E meix , 
though one has to keep in mind that all limitations listed 
above. As long as E max = O(10) TeV (within a factor a 
few), its exact value is of minor impact for predictions of 
p/p at E < 1 TeV. Qualitatively, a higher value of E max 
would increase the slope of the rise in the positron ratio, 
while a lower value would flatten it. 

Finally, we comment on the role of nuclei in our calcu- 
lations: in [l5| it was found that for a typical composition 
mixture like the one measured locally, correcting for this 



effect roughly amounts to a factor e ~ 1.20. We have 
repeated the effective weighted-average renormalization 
using the same cross-section weights as in [l5j], but the 
updated composition ratio compiled in [lfjj], table 24.1. 
We obtain a factor e =1.26, which we shall use to renor- 
malize both Eq. © and the ISM contribution. The ISM 
spallation contribution to p/p can be written in the same 
approximation as above in the form 



J p jsm{E) 2eX(E) 



m v E^-i- & JE 



J P (E) 

with the grammage parameterized as: 
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JE 



(8) 



x(E) = r 
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I gem" 2 (£>10GeV). (9) 



Here we adopt 5 = 0.6 and T — 5.5, well within the range 
discussed in UM- 
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FIG. 1: The p/p ratio for the parameters reported in the 
text, together with a simple model of secondary production 
in the ISM (dashed line), and with the recent data from 
PAMELA [llj]. The dotted and dot-dashed lines represent 
the contributions of the A and B terms of Eq. (f2J alone, re- 
spectively. The thick solid curve is the overall contribution 
due to ISM plus the new mechanism, while the thin solid 
curve only includes the ISM contribution plus the B term. 

The predictions thus obtained are reported in Fig. [Tj 
Although the model described is very simple, the overall 
agreement with the data is good, with the predictions for 
the conventional model (only antiprotons from spallation 
in the ISM) strongly differing from the present ones be- 
yond the E ^100 GeV region. In the case considered 
here the ratio flattens at first, then grows with energy. 
The latter behavior is exclusively due to the .A-term of 
Eq. while the former behavior is due to the interplay 
of the decreasing ISM term of Eq. ([8|), the rising .A-term 
and the relatively flat B-term of Eq. ([5]). It is impor- 
tant to note that the S-term accounts for production of 
p without "reacceleration" , thus it does not depend on 
the diffusion properties, only on the density of the envi- 
ronment ni and the typical timescale t$n'- we see that 
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its presence alone contributes to change appreciably the 
shape of p/p at E >100 GeV. The role of the ,4-term 
is even more dramatic at high-energy (and indeed it is 
essential to explain a positron fraction rise of the kind 
revealed by PAMELA), but it is somewhat more model 
dependent both in shape and normalization. 

Discussion and Conclusions — In this letter, we have 
discussed an important signature of the mechanism pro- 
posed in [l[ to explain the anomalous behavior of the 
positron ratio at high energies: a harder component 
should emerge in the antiproton spectrum at energies 
above ~ 100 GeV. The p/p ratio flattens at first, then 
eventually starts rising with energy. New data at high 
energy from PAMELA and AMS-02 [ig[ should easily 
distinguish between this explanation and a pulsar related 
one for the positron fraction. As discussed above, though 
the effect predicted here (and in [l| for positrons) must 
be present, its strength depends on the many parameters 
of the problem and on whether they are appropriate to 
describe the final stages of SNRs. This uncertainty is 
mainly of concern for the rising ( "reacceleration" ) term, 
while the injection term is less model dependent. The 
latter shows-up as a flattening in the antiproton ratio 
and represents a conservative prediction for the energies 
just above the ones currently probed by PAMELA. Even 
limiting ourselves to the effects of the latter term, the 
implications for astrophysics are of crucial importance: 



The good news is that the high-energy range of the an- 
tiproton spectrum may reveal important constraints on 
the physics of the CR acceleration sites. The bad news 
is that it is not straightforward to infer from high energy 
p/p-data the propagation parameters, as the diffusion in- 
dex 5, since they are partially degenerate with source pa- 
rameters: The thin solid line in Fig. [TJ might be easily 
confused with a purely ISM model with no contribution 
at all from SNRs, but a lower value of the diffusion index 
in the ISM, 5. 

Similarly, our results may change dramatically the per- 
spectives for the detection of DM via a signature in high- 
energy antiprotons: Indeed, we have discussed a purely 
astrophysical mechanism to produce a high-energy "ex- 
cess" of antiprotons over the secondary yield from ISM 
production. Even a subleading role for the mechanism 
proposed in [l[ in explaining the positron excess might 
produce measurable anomalies in the antiproton spec- 
trum. An "excess" in the high-energy range of p/p could 
not be interpreted anymore uniquely as manifestation of 
new physics: compare for example Fig. 3 in [l2| with 
our Fig.[TJ The mechanism proposed here might thus re- 
quire a paradigm change for DM searches via anti-matter, 
at least until the contribution from standard astrophysi- 
cal sources is understood and corresponding uncertainties 
are kept under control. 
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